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We have recently proposed an efficient computation method for the frictionless linear elastic 
axisymmetric contact of coated bodies [A. Perriot and E. Barthel, J. Mat. Res. 19 (2004) 600]. 
Here we give a brief description of the approach. We also discuss implications of the results for 
\Q ' the instrumented indentation data analysis of coated materials. Emphasis is laid on incompressible 

or nearly incompressible materials (Poisson ratio v > 0.4): we show that the contact stiffness 
rises much more steeply with contact radius than for more compressible materials and significant 
elastic pile-up is evidenced. In addition the dependence of the penetration upon contact radius 
O , increasingly deviates from the homogeneous reference case when the Poisson ratio increases. As a 

result, this algorithm may be helpful in instrumented indentation data analysis on soft and nearly 
incompressible layers. 
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I. INTRODUCTION 



-a 

Understanding the elastic contact to a coated substrate is a prerequisite to the analysis of more complex thin film 
• behaviour. It is also an interesting problem because of two difficulties: 

1. the coupling by the elastic field of planar parallel interfaces, which determines the response function; 

2. the mixed boundary conditions usually involved in contact problems. 

The problem has been dealt with in numerous publications, for example 0, 13, 0, H| ■ In a recent contribution, we 
have proposed a numerically efficient approach [||. We first summarize this technique, insisting on the structure of 
lO ' the elastic response function on the one hand and on the treatment of the mixed boundary conditions on the other 
hand. 

We also apply the resulting algorithm to the specific case of nearly incompressible or incompressible layer material. 
Compared to compressible layer material, the stiffness increases much more rapidly. At the same time, the dependance 
of the penetration upon contact radius is also strongly affected. More generally, the deviation of the penetration vs. 
contact radius relation from the homogeneous substrate solution could limit the direct applicability of an Oliver-Pharr 
approach to instrumented nanoindentation data treatment for coated substrates. 



I 

O . II. CONTACT TO A COATED SUBSTRATE 

O 



The present method is based on the elastic response function of a coated substrate calculated by Li and Chou 0. 
The expression of the response function is complex but is characteristic of coupled planar parallel interfaces. We will 
first highlight these main features. 

A. Response functions of coupled interfaces 

The coupling of parallel planar interfaces through the layer of material 1 (thickness t) is schematized on Figure ^ 
In-plane symmetry suggests Fourier transform in the x,y plane (wave vector k). If the fields obey an equation of the 
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laplacian type, then the response function D obeys, in the notation of Fig. ^ 



d 2 

w 2 - — ) D(k, z) = for z < -t (1) 
d 2 

w\ — r—T ) D(k, z) = A(k)S(z - z') for - t < z < (2) 
dz z 

d 2 

dz 



w l ~ — ) D ( k , z ) = for < z (3) 



with adequate boundary conditions on the interfaces and at z = ±oo. A(k) is an excitation located on the plane 
z = z' . For a harmonic field n = 1, and n = 2 for a biharmonic field. This general structure will be specified below 
along with the parameters Wi. 



1. Electromagnetic coupling between planar interfaces 



The interaction of the electromagnetic field radiated by polarisable bodies at close distance results in the van der 
Waals interaction. Assuming that two such media with dielectric constants and parallel interfaces are separated by 
a distance t, the field fluctuations [S| with wave vector k in the x, y plane and angular velocity uj n obey Eas lllHl with 



Wi = y/e.col + k 2 (4) 

The response function can be calculated. For instance, the response at z = to an excitation with wave vector k 
located at z — (local response) is 

D{0,0)—~ (5) 

with 

V = 1 - exp(-2wii)A iAi2 (6) 

and 

A. = ^ (7) 

This form of the response function is characteristic of coupled parallel interfaces. Indeed the denominator T> is the 
determinant resulting from the solution of the linear system Eqs. 1 1131 It couples the exponential field propagation 
factor in medium 1 (i.e. exp(— 2w\t)) with the polarisability mismatch coefficients Ay . These antisymmetric mismatch 
coefficients result from the boundary conditions at the interfaces, and appear quite generally in interfacial problems 
such as electrostatic image potential, Fresnel reflection coefficients or Dundurs elastic mismatch parameters. 

The response function D can be used to calculate van der Waals interactions. Indeed the secular equation T> = 
determines the coupled surface plasmons eigenfrequency from which the variation of the energy of the system with 
the distance t can be calculated. 



2. Response function of a coated substrate 



Similarly for our contact problem, we assume a coated elastic half-space under an axisymmetric frictionless loading 
(Fig. |2J). The layer, of thickness t, and the half-space are elastic, isotropic and homogeneous, while their adhesion is 
supposed to be perfect. Let Eq and Vq (resp. E\ and v\) be the elastic modulus and the Poisson ratio of the half-space 
(resp. of the layer). 

Due to axisymmetry, we use / the O^-order Hankel transform of / defined as: 

/•oo 

7(fc) = / drrJ Q (kr)f(r) 
Jo 

instead of the Fourier transform. Here Jq(x) is the O th -order Bessel function of the first kind which displays a 
cosine-like oscillatory behaviour suitable for the present geometry. 
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Then, for a coated substrate, we fall under the framework of Eqs. 11131 for linear elasticity. As is usual in 2-D 
elasticity problems, the fields obey a bilaplacian equation so that n — 2. We consider a static problem so that Wi = k. 
Under these hypotheses, Li and Chou obtained a relation Q between the applied normal surface stress q(k) (taken 
positive when compressive) and the displacement u z (k,z) (positive is inwards). 

More specifically for z — (with u z (k,0) = u(k)): 



ku(k) = C(kt)q(k) (8) 



where the response function is: 



an) = a i+4bfcte ' 2fc ;- abe - 4fct (9) 

E{ 1 - ( a + b + Ab{ktf)e- 2kt + ab e" 4fct 

-71 h — n — - E l( 1 + t/ o) 



b= ^7, <*=%£Bk,1i =3-4^ and 70 = 3 -4^ 



l+a-fo ' " a+7l' Eoil+vt) 

El being the reduced modulus of the layer defined as E\/{1 — v\ 2 ). 

The structure of the response function is again characterized by the exponential field propagation factors and anti- 
symmetric elastic mismatch factors for the field transmission at the interfaces. However, compared to electromagnetic 
response, the form is complexified by the bilaplacian-type field and the shear-free boundary conditions at the surface. 



B. Boundary conditions 

Contact problems such as indentation are characterized by mixed boundary conditions (Fig. [2J : the normal surface 
displacement u is specified in the contact zone - by the penetration and shape of the indenter - while the normal 
applied stress q is given outside - usually zero -. Now Eq. JSJ can only be used to calculate the surface displacement 
everywhere on the surface provided the normal surface stress is known everywhere on the surface. The response 
function Eq. [21 is therefore useless as such pfl- However, this problem can be circumvented by a mathematical trick 
(or change of basis, or use of auxiliary functions). 

The method, which appears explicitly in the later papers by Sneddon 0, Q , has proved fruitful for the description 
of complex contacts. For example we have been able to propose an exact solution to the long standing question of 
the adhesive contact of viscoelastic spheres in this way [lOj . 

Algebraically the method relies on the relation 

Jo Vf - 1 

that is to say the Bessel function Jo is the cosine transform of the function 

9(r-t) 



Vr 2 - t 2 



(11) 



where 0, the Heaviside step function, anticipates our mixed boundary conditions. 

The trick can also be viewed as transforming back by a cosine transform the variables which were transformed 
forward with a Hankel transform. We do not fall back in proper real space, but in a space where our boundary 
conditions assume a suitable form, due to Eq. ^3 

Explicitly we introduce the auxiliary fields g and 6 defined as ^(| : 

/>oo 

g(s) = / dkq(k) cos(fcs) (12) 
Jo 



9(s)= dk ku(k) cos(/cs) (13) 



Let us now apply the cosine transform to Eq. (|SJ). Then: 



e(s) = -l g(r)[l dk C(kt) cos(kr) cos(ks) I dr (14) 



7T 
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Similarly, Eq. El becomes 



f +co sq(s)ds 

9{t) = I 7W=W (15) 



As a result of the boundary condition q(s) = for s > a, g(t) = for t > a where a is the contact radius. Then the 
upper bound in the spatial integral Eg. 1141 is actually a, not infinity. The equilibrium equation 1141 has now the form 



0{S) = lf 
* Jo 



g(r)K(r,s)dr (16) 



and a simple algorithm can then be devised. 

Indeed the function 9 is known inside the contact zone for an arbitrary indenter shape h(r): from Eos. 1101 and 1131 
one obtains (s < a) 

„, ^ d f s , rh(r) 

" */„ *7fH» (17) 

This can be calculated analytically for most simple indenter shapes. In addition, the kernel K in Eq. El a s defined 
by Eq. El can be calculated from the response function C(kt) by fast Fourier transform (FFT). 

Thus, discretizing Eq. 1141 g(s) for s < a can be calculated numerically as the solution to a linear system. From 
g(s), s < a, force, penetration and contact stiffness can be computed [||. 

Some of the results obtained with this algorithm have been presented earlier. In this paper, we want to insist on a 
specific case: when the layer is incompressible or nearly incompressible. 



III. INCOMPRESSIBLE AND NEARLY INCOMPRESSIBLE COATINGS - EFFECTIVE CONTACT 

STIFFNESS 

In our previous paper |(J, we calculated the contact force P, the contact stiffness S and the penetration S as a 
function of contact radius a for a large range of modulus mismatch and a single value of Poisson ratio v = 0.25 for 
both substrate and film. These results provided a detailed description of the transition between the film dominated 
regime and the substrate dominated regime. 

It is well known, however, that confinement of the layer at large contact radius values a/t 3> 1 will result in specific 
phenomena for nearly incompressible materials ^lj • Indeed, for such materials, volumetric deformations are penalized 
so that shear deformations predominate. In the present case, such deformations are hampered by the confinement. 
The response is therefore dependent upon the axial compression modulus in the absence of lateral strain - the so-called 
cedometric modulus of the layer - 

E{l- V ) 

E ° - (l-2„)(l + „) (18) 

Especially relevant is the case of a compliant layer deposited on a more rigid substrate. This could be a polymer 
film on a glass substrate for instance. Using our previous normalization scheme 0, we have calculated the response 
of the elastic frictionless contact of a sphere on a coated substrate for a film/substrate reduced modulus mismatch 
equal to 0.1. The numerics were carried out with the Igor data treatment software (Wavemetrics) on a standard 
1.60 GHz processor. To probe the relevant part of the response function, a variable cut-off depth B = 3000a /t was 
used. The other parameters for the numerics are Ni = 2 20 points for the cosine transform and the contact radius 
a was discretized over N2 = 500 points. The resulting computation time for a given contact radius is of the order 
of 1 second. The results at large Poisson ratios were checked by increasing the computation parameters Ni and N2 
without significant variations in the results. 

Fig. |3] displays the normalised contact stiffness, or effective reduced modulus, as a function of the normalized 
contact radius for Poisson ratios comprised between 0.1 and 0.5. The usual transition between film modulus (0.1) and 
substrate modulus, which is the normalizing parameter, i.e. 1, is readily obtained. The results are almost insensitive 
to the Poisson ratio v for v < 0.25. Sizeable deviation is recorded for 0.4 < v < 0.5. In this range, we observe that 
the transition occurs earlier and is increasingly steep as the layer becomes less compressible. However, the system 
does finally reach the bare substrate reduced modulus. 

For incompressible coating materials, this earlier saturation of the reduced modulus is essentially due to the elastic 
compliance of the substrate itself and is better evidenced from the surface and interface normal displacements. These 
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displacements as computed by inverse Hankel transform from the discrete g functions obtained at various a/t ratios 
are displayed on Fig.0] At low penetrations (a/t = 0.1), the substrate remains undeformed and the compressible and 
incompressible systems are undistinguishable. At larger contact radii (a/t = 1), the substrate starts to deform and the 
coating material is increasingly confined. For the incompressible material, an elastic pile-up of the coating material 
forms around the contact zone, while a global elastic sink-in results from the deformation of the substrate and increases 
with increasing contact radius. Generally speaking, these results illustrate the impact of the incompressibility on the 
partitioning of the elastic displacements between substrate and coating. Most noteworthy, at contact radius a/t = 8, 
for an incompressible layer, the surface displacement is almost completely due to substrate deformation; yet a roughly 
equal contribution from layer and substrate is calculated for a compressible coating material under the same conditions 
of confinement. 



IV. SURFACE DEFLEXION AND INSTRUMENTED INDENTATION DATA ANALYSIS 

With instrumented indentation, in the absence of direct measurement of the contact area, the Oliver & Pharr 
method is used to infer the contact radius from measured variables, i.e. force, penetration and contact stiffness |l3j |. 
One usually uses 

hc = 6- e| (19) 

where S is the contact stiffness, P the force, 8 the penetration and h c the contact depth. The form of the equation 
and the value of the constant e are such that eq. ^5] is exact for a purely elastic system. It has been experimentally 
demonstrated that eq. EH is very useful also for elasto-plastic materials, mostly when plastic pile- up is minimal. 

For the analysis of instrumented indentation on coated substrates, the same equation is used. However, Ea. 1191 is 
no longer valid. Similarly, the usual relation between penetration and contact radius such as 

5=1— a (20) 
2 tan/3 v ; 

for a cone of half-included angle breaks down. It is replaced by a more complex relation which involves the mechanical 
parameters of the system. This explicits the fact that Eq. 1201 is a mechanical relation even though the mechanical 
parameters have accidentally dropped out for a homogeneous substrate. For coated substrates, some examples are 
plotted on Fig. [5j A salient feature is that small contact radii are more easily reached ( i.e. with smaller penetrations) 
for incompressible layers than compressible layers. This is due to the elastic pile-up effect. However, for a/t ~ 4.5 
a cross-over takes place. Larger radii are achieved by larger penetrations for incompressible layers because of the 
increased effective stiffness. 

If we keep the form of Ea. 1191 for coated substrates, we may attempt a more accurate approach by calculating the 
value of e for a given coating configuration. Some results are displayed on Fig. Note that in the transition region, e 
drops to significantly lower values and that again the approach to incompressibility results in considerable variations. 
Values for e as low as 0.2 may affect the data treatment in cases where the elastic contribution is significant (large 
S, ie sizeable elastic recovery). Due to the ease of calculation involved here, adequate e values are readily calculated 
self-consistently during actual data treatment. 

Finally, it is important to observe that failure to reach a correct evaluation of the contact radius is a double source 
of errors. A direct error incurs on the evaluation of the effective reduced modulus E* through 

E* = (21) 

where (3 is the non-axisymmetry correction factor and A the contact area. A second error affects the evaluation of 
the layer reduced modulus from the effective reduced modulus because the model for the latter, such as discussed in 
section ITTT1 will be evaluated at an incorrect value of the contact radius. In the case of compliant films, these errors 
add up. An evaluation for v = 0.2 and a modulus mismatch of 100 results in an error of up to 60 % at a/t ~ 3. 



V. CONCLUSION 

Summarizing our approach to the elastic contact of coated substrates [(j, we have shown that the form of the 
response function calculated by Li and Chou |jj is characteristic of the coupling of parallel planar interfaces through 
the elastic field. Using a new basis of functions well suited to the mixed boundary conditions, we obtain an integral 
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relation that is easily solved numerically and provides force, penetration and contact stiffness as a function of the 
contact radius, for arbitrary moduli ratios and arbitrary axisymmetric indenter shape. 

The stiffness of incompressible layers increases much more steeply with contact radius than more compressible 
coatings as soon as the layer is confined in the contact. Indeed volumetric deformation is strongly penalized while 
significant shear is prevented by the confined geometry. 

This is accompanied by sizeable elastic pile-up which significantly reduce the penetration needed to reach a given 
contact radius. For large contact radii, the stiffness of the layer dominates and the surface deformation is almost 
completely due to the elastic yielding of the substrate. Again, this phenomenon occurs at lower contact radii for 
nearly incompressible layers because of the enhanced stiffness. 

The deviation from the homogeneous material penetration vs. contact radius relation presumably bears upon the 
validity of the Oliver and Pharr approach to instrumented indentation. A possible improvement for soft and quite 
elastic layers could be obtained by calculating more accurate values of e with the present method. 
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CAPTIONS 

Fig. 1: Geometry of the coupling between parallel planar interfaces. The layer thickness is t. 
Fig. 2: Schematic representation of the indentation of a coated elastic half-space. 

Fig. 3: Effective reduced modulus S/2a as a function of normalized contact radius a/t for a layer/substrate reduced 
modulus mismatch of 0.1 (E* ubstrate — 1, Ej ilm = 0.1) and different layer Poisson ratios. The substrate Poisson 
ration is 0.2. 

Fig. 4: Computed normal displacement u of surface (plain) and interface (dashed) for the contact of a sphere to 
a coated substrate (E* ubstrate / EJ ilm = 10). Radii are normalized to contact radius a, displacements to a 2 / R. The 
coating Poisson ratio is v = 0.2 on the left, 0.5 on the right. Calculation parameters as in Fig. |21 The curves are 
offset by unit increments for clarity. 

Fig. 5: Computed penetration as a function of normalized contact radius for a sphere on a coated substrate 
(E substrate / Efiim = 10' ^substrate = 0.2) for various film Poisson ratio. The penetration is normalized to a 2 /R. 

Fig. 6: Computed values for the "constant" e in Eq. EH as a function of normalized contact radius and film Poisson 
ratio. Same parameters as in Fig. 5. 
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